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16 . Abstract 

Carried out herein is generalization of the flow pass 
method for a stationary problem of diffusion of a minor 
component in a stratified medium, with the presence of the 
force of gravity and. chemical reactions. 

It is shown that the modified flow pass, method, applied to 
problems of chemistry of planet atmospheres, possesses con- 
siderable effectiveness, both in the case when the coeffici- 
ent of diffusion changes severely in the examined region, 
and in the case when diffusion is the prevalent process, . 
as. compared with chemical reactions, i.e . . in. the case when 
a regular pass proves inapplicable, or applicable in a 
limited interval of the decisive parameters. 
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MODIFICATION OF THE FLOW PASS METHOD AS APPLIED TO PROBLEMS 
OF CHEMISTRY OF PLANET ATMOSPHERES 

V. A. Parshev 


Carried out herein is generalization of the flew pass method 
for a stationary problem of diffusion of a minor component in a 
stratified medium, with the presence of the force of gravity and 
chemical reactions. 


It is shown that the modified flow pass method, applied to 
problems of chemistry of planet atmospheres, possesses considerable 
effectiveness, both in the case when the coefficient of diffusion 
changes severely in the examined region, and in the case when 
diffusion is the prevalent process, as compared with chemical 
reactions, i.e. , in the case when a regular pass proves in- 
applicable or applicable in a limited interval of the decisive 
parameters . 


The problem of determining the spatial distribution of 
minor components in the atmospheres of planets amounts to the 
solution of a system of equations, each of which has the form 



with one of the following typical boundary conditions at the 
upper boundary: 

al / — representation of flow, 

*2/ |z.s2Af * — repi asentation of concentration 


and at tne lower boundary: 


♦Numbers in the margin indicate pagination in the foreign text. 


— representation of flow, 


“'"•L-frAt 

»3/ |z »0 Sa ttiO 


representation of photo 
chemical equilibrium 


— representation of concentration 


Here, n i is the concentration, and 4^ is the flow of the 
component i. 




P i( L i are the rate of production and the function of particle 
destruction, k is the coefficient of diffusion, brought about 
by motion of the medium, ^ is the coefficient of molecular 
diffusion, H if H cp are the scales of the altitudes of the examined 
component and the medium, respectively, T is the temperature, 

Oj s the thermodiffusion factor, Z is the altitude reckoned from 
a c. 'tain level (not necessarily from the surface of the planet, 

0" ,ZM) , ZM is the upper boundary, and i is the number of the 
component . 

Using the method of component breakdown, the problem of 
many components can be reduced to the sum of the problems of 
the given type, which makes it possible to examine the problem 
(1.1) -(1.4) for a single component. Subsequently, the index i 
will be omitted. 

The following boundary value problem corresponds to the 
differential boundary value problem (1.1) -(1.4)* 


2 


- ai*uW«.-< + a 1 . ( + £;)*; -g. n .. <=dii ( Z .i) 

Fi = - ( a; H - tf t v, « c ) (tiUH-i) (It) 
ai/ />i» = Qa i > a 2/ rt-iu = *1 a, , (Z.5) 

il/ F,= Q e ; H/ n,= d,/£ s ; id/ n< = n 0 . (&4j 


The equation is approximated according to three points on a 
uniform network, according to Z, in a conservative manner, u is 
the number of the point in sampling for Z, ^ is the difference 
analog of L, d A is the difference analog of p, and F i is the 
difference analog of the flow <j>. 

By bounding our problem by the description of the method of 
solution of the difference system, according to all questions 
associated with the derivation and substantiation of the differ- 
ential equations, as well as the switch to an approximating 
system of difference equations, we will refer the reader to the 
appropriate literature (for example, one can find a detailed 
bibliography on these questions in study [5]) • 

In order to solve problems of the type (2), the pass method [13 
is usually utilized, which possesses great effectiveness and 
stability and makes it possible to obtain a solution of the problem 
in a rather broad range of decisive parameters. However, there 
exists a large class of problems of the described type, when the 
use of the standard pass method (SP) is not possible, since it 
leads to a considerable loss of accuracy, for example, in problems 
of chemistry of planet atmospheres: (*) when the coefficient of 
diffusion changes severely, i.e. , the coefficients a. } , b in (2) 
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charge in the examined area by several orders} or (**) , when 

the diffusion is a more effective process, as compared with the 

chemical reactions between the components, and the flow through 

the boundary is absent, i.e. , a i l( b^ and Fj.=I^ z =0. Under 

these conditions, the total concentration of the component in 

the examined area is determined by the balance between the total 

production and destruction, and its altitude distribution is 

brought about by diffusion. The use of the standard pass method 

for solving these problems leads to a considerable loss of 

accuracy, if the order of the relationship of the minimum value 

of b » e , d, , a. to the maximum value of a. ,,b. is greater 
1111-I 1-I1 

than the number of significant digits with which the computer 
operates. The flow pass method (PP) [2,3,4]] makes it possible 
to obtain a solution of the problem in both cases (*) and (**) , 
if the flow is determined as 4> =k-^B. • 


Given in the present study is a modification of the flow 
pass method for a flow of a more general type (1.4), which makes 
it possible to overcome both of the indicated difficulties. 


Description of the Algorithm 

Thus, we will examine the difference boundary value problem 
(2.1)-(2.4). For a common pass, we have 


+ , = 2 ; n t +m. ; , 


Z.-< = 


j i 



f 



m-i-, = 


di ± m.- a / 


( 3 . 1 ) 
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The first equation (2), we will write in the form 
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^ “■ *" * ^4 "t" di * 

Expressing n. i =-•£-*• + i -**n. from (2.2) and substituting 

i + ^ a i 1 

into (3-1) » we obtain 

-Fi+n.1 (£>, - Zi&i) =: n^a.;. (5) 


Introducing new coefficients according to the formulas 


9 °^i = fO 

we will write (5) in the form 

- Fi = - dint 1-&V (?) 

From (?) and (4), excluding n , we obtain the recurrent relation- 
ship for the flow 



oil 




z-^JL 
+ 6 ^ 


From (6), we obtain the expressions for 2i , m 1 

H- = ■ w ,= 


From the first equation in ( 3) * we obtain the recurrent relation- 
ship for n ji 










The recurrent relationships for 0t i *Y i are obtained from (6), (3*2), 
and (3-3) * 
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( 3 ) 



gLc +4 ^ 

<X t + a ; 


t< 


With regard for (9)* expression (8) takes on the form 

ft;*, = (C< «•; + <v,-*-K„)/ + +<Z 0 - 


The final formulas for the direct pass have the form 



oC 1 

. y — ct >t (jq A 

cLi- 

&£+«».,+ £i« 

> 0 i a t + oL*~<+£<*< * 1 • ' 

a </ 

= 0 > ¥nz 

= imi 


di ,v ~ ^ni 



/ 


Reverse pass is carried out according to the formulas 


Hi-M - 71. ' i7^^7, n ‘ + 


d,„ + 


c — - — /C. + f'w^vrf.v. y.v, 

6, », ^ ‘*•♦1 * <^i*i 

f, 5o > vii = (jfi+Qo) /oti 9 (ii) 

(i/ fti fa > 

£ 3/ /J { ' 7 £/ - °*1 ' “ ^4 . 


The obtained formulas, as easily noted, are similar to the 
formulas of flow pass, given in [ll? therefore, they may be 
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called the algorithm of the modified flow pass (MPP) . 


The stability of the recurrent relationships given above is 
evident. 


Discussion of the Method 

We will give reasons which show that the calculations 
according to the given scheme prove more accurate than those 
which are carried out by the traditional pass method. 


Insofar as the coeff icientsaj, b ^ ej* and d A are positive, 
the pass coefficients a i* Y i > 0, and, thus, the sums of the 
positive magnitudes stand both in the numerator and in the 
denominator of the expressions (10) , and, although the absolute 
error in the sum is made up of the absolute errors of the terms, 
the relative error does not increase, insofar as the moduli of 
the magnitudes in the numerator and denominator increase, while 
the accumulation of the error (relative) during division and 
multiplication does not have such catastrophic consequences as 
during subtraction. 


If d. , e ,«a. , d and a =0j y =0, then the method of mathe- 

111-I1 1 z nz • 

matical induction may be used to show that a and y , 41 d 

n i i i*l t.j k ’ 

i .e . , magnitudes of a single order take part in the determination 
of a i , » although only as sums, and the relative error does 

not increase. 
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These reasons show that the modified flow pass is free from 
the above-indicated difficulties, which occur during the solution 
of problems of chemical kinetics by normal pass. 

For illustration, we will examine two models of the problem. 
The first is described by the following system of equations 



d $ _ Q 

- 2> 

(b) 

=1 > 

2) = exp ); 


“t 5 ! ZA. “ C , 

ze iC : 
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where 3>0 is the parameter of the problem. 


There exists an accurate solution of the prob?.em, which has 
the form 

n. <?* p (-2). 


The second problem is formulated using the following system 
of equations 



■V# 

Ct 2 


E (i-n) , 







where the coefficient of diffusion is assumed as equal to 1, 

E« 1, ZM=20. The asymptotic solution of this problem, with the 
condition ZM*E «1, has the form 

ft “ ft. c p 2 j , where — 2 M . 

For these problems, in accordance with what has been set 
forth above, approximating difference equations were written, 
which were solved both by the method of the standard pass and 
by the method of the modified flow pass. A computer, operating 
with 12-digit numbers, was utilized for the solution. 

Problem (A) was solved with various values of the parameter /g 
3 in the coefficient of diffusion in a uniform network with a 
spacing a Z= 0 .1 . It turned out that the method of the standard 
pass, with 3<2, provides a solution of the problem which differs 
by less than 2# from the analytical solution, but with 3>2, the 
solution differs considerably from the analytical solution (fig. 2). 

The method of the modified flow pass provides a solution which 
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coincides, with an accuracy of up to 2$ for all 0 (fig. 1). 

Problem (B) was solved with various values of E in the very 
same network. It was found that, with E<10“ 10 , the solutions for 
the standard pass and the modified flow pass coincide, with an 
accuracy of up to 5$. but with E>10“ 7 , these solutions differ 
from the asymptotic solution (fig. 3)# and only with E<10“ 7 do 
all three solutions coincide — that obtained according to the 
standard pass up to 10$, and that obtained according to the 
modified flow pass with an accuracy of up to 4$ (fig. 4). 

With E<10~ 10 , the solution by the method of the standard 
pass is impossible to find (the system becomes confluent for 
the computer) , while the modified flow pass provides a solution 
with extremely small E, for example, a solution was obtained with 
E=10” ; 0 (fig. 4). 

The results given above show the effectiveness of the method 
of the modified flow pass during the solution of problems of the 
described type, for example, those associated with the chemistry 
of planet atmospheres. It combines within itself both the ad- 
vantages of the standard pass, namely efficiency and a small 
volume of required storage, and the advantage of methods which 
utilize a double length of the computer number — namely, high 
accuracy of solution. 

In conclusion, the author th?nks V . A. Krasnopol’skiy for 
his detailed discussion of the study. 
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Fig. 1. Dependence of relationship^ (n/r^ 0 on the 

coordinate Z. (T) n=r^ na | (2) n is the solution 
of the problem Dy the method of the modified flow 
pass and standard pass with S<2. 


Fig. 2. Dependence of magnitude of ^(n/r^naj ) on the 
coordinate Z for problem (A) by the method i 

Modified flow pass for all Pj (2) Standard pass 
with P = 2.5j(jf) Standard pass with P=3. 
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Fig. 3- Dependence of magnitude of £n(n/n aS y mp ) for 
problem (B) on the coordinate Zt 

(l) E=10" 5 ; solutions by modified flow pass and 
standard pass methods coincide; 

(D Asymptotic solution and solutions by methods of 
modified flow pass and standard pass with E=10“ 1 0 
coincide . 



Fig. 4. Dependence of magnitude of (n/n asymp ) on coordi- 
nates Z lor problem (B):(I> Asymptotic; (2) Modified 
flow pass with 10~ 8 jE ^lO" ' °; (3) Standard pass with 
E=10~ 10 . 


